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D ■ Under pressure cesium undergoes a transition from a high-pressure fee phase (Cs-II) 

-J. , to a collapsed fee phase (Cs-III) near 4.2GPa. At 4.4GPa there follows a transition to 

the tetragonal Cs-IV phase. In order to investigate the lattice vibrations in the fee phase 
and seek a possible dynamical instability of the lattice, the phonon spectra of fcc-Cs at 
volumes near the III to IV transition are calculated using Savrasov's density functional 
linear-response LMTO method. Compared with quasiharmonic model calculations including 
non-central interatomic forces up to second neighbours, at the volume V/Vq = 0.44 (Vq is 
the experimental volume of bcc-Cs with ao=6.048A), the linear-response calculations show 

'^ ! soft intermediate wavelength Triigi [^^0] phonons. Similar softening is also observed for short 

a ■ 

O ' wavelength -L[^^^] and i^[00^] phonons and intermediate wavelength L[^^^] phonons. The 

Born- von Karman analysis of dispersion curves indicates that the interplanar force constants 

exhibit oscillating behaviours against plane spacing n and the large softening of intermediate 

wavelength T^ijo] ['C'CO] phonons results from a negative (llO)-interplanar force-constant $ri=2- 

The calculated frequencies for high-symmetry K and W phonons and longitudinal X and L 



C/3 



C^ 



CJ 
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Q ■ phonons decrease with volume compression. In particular, the frequencies of the T'[iio] [^^0] 



O 
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X 



phonons with ^ around | become imaginary and the fee structure becomes dynamically 
unstable for volumes below 0.41Vo. It is suggested that superstructures corresponding to the 



r^ , qT^O soft mode should be present as a precursor of tetragonal Cs-IV structure. 
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1 Introduction 

At ambient pressure cesium metal crystallises in a bcc structure with experimental lattice 
parameter a^ = 6.048ApP at OK, obtained from the high-temperature values by extrapo- 
lation. It is a simple s^ metal with nearly free-electron character. Under pressure the s 
valence electrons are transferred to more localised d-like statesj^, and cesium exhibits an 
interesting sequence of phase transitions [3-8]. At the pressure 2.3GPa a transition [HlllI from 
bcc (Cs-I) to fee (Cs-II) occurs with a small reduction of volume. Near 4.2GPa, the fee Cs-II 
phase undergoes an isostructural transition ^j to a collapsed fee phase Cs-III with a large 
volume reduction (9%). Then at 4.4GPa there follows a transition to the tetragonal Cs-IV 
phase [S] with a decrease of the coordination number from 12 in Cs-III to 8 in Cs-IV (strictly, 
4 nearest neighbours at 3.349Aand 4 second- nearest neighbours at 3.542Aat 8 GPa). At 
higher pressures more complicated phases [HI 13 E] appear. 

Many detailed investigations have been devoted to the study of the phase transitions 
in cesium, both experimentally (see Ref.jZj and literature cited therein) and theoretically 
(see, e.g., Refs. [7-12]). The pressure-induced electronic s-^d transition was believed to be 
the driving force for destabilising the highly symmetric low-pressure structures (bcc or fee) 
with respect to lower symmetry structures. The transition III— i>IV was attributed to the 
accelerated progress|3] of the s^d transition and the unusual decrease of the coordination 
number from Cs-III to Cs-IV has been interpreted in terms of directional bonding induced 
by the s^d transition|12| IT5] . 

In general, the s^>-d electronic transition in cesium affects not only the static lattice prop- 
erties by changing the electronic band structure but also the lattice vibrational properties 
by modifying the effective interatomic interactions. The structural behaviour of cesium may 
therefore be reflected in the phonon spectrum and it is therefore of interest to investigate 
the phononic behaviour with applied pressure. Glotzel and McMahan^Tj suggested that an 
anomaly in the lattice vibrational contribution to the pressure is a possible mechanism for 
the isostructural II— i>III transition. Recently, Christensen et al.[16| calculated phonon dis- 
persions of fcc-Cs within the quasiharmonic approximation to study the thermal expansion 
coefficient of cesium and the isostructural transition. They found that below V/Vq = 0.375 
(Vq is the experimental volume of bcc-Cs at ambient conditions) fcc-Cs becomes unstable 
due to softening of a transverse [110] phonon branch, which is characterised by a negative 
shear elastic constant C. Very recently, Xie et al. J7] investigated the phonon instabilities 
of cesium using density-functional perturbation theory. Their calculations showed that the 
instability of the fee cesium occurs at a volume of about 0.46Vo, which is also driven by a 



negative C . However, it is not clear what happened with the phonons before the C goes 
soft. Using the density functional linear-response LMTO method jTHj, we calculate the accu- 
rate phonon spectrum of fcc-Cs at the volumes near the III-IV transition to study the lattice 
vibrations in fcc-Cs and to seek a possible dynamical instability of fcc-Cs before C becomes 
negative. According to our calculations, fcc-Cs becomes unstable at a volume slightly smaller 
than V/Vq = 0.41. Soft Tjijoj-f^^O] phonons with ^ ~ | are found to be responsible for the 
dynamical instability. 

This paper is organised as follows. Some computational details concerning linear-response 
LMTO calculations are described in the next section. Section 3 presents the results ob- 
tained and some discussions. Firstly, the phonon dispersion curves for fcc-Cs at the volume 
V/Vq = 0.44 are analysed in terms of force constants between atoms; then the phonon fre- 
quencies at various volumes are discussed in relation to the instability of the fee structure. 
Several concluding remarks are given in the last section. 

2 Computational Details 

The density functional linear- response LMTO method [18j is used in the present study to cal- 
culate the phonon spectrum of fcc-Cs at the volumes around the Cs-III - Cs-IV transition. 
Phonon dispersion curves for a large number of simple and transition metals and a few com- 
pounds have been calculated by the linear-response LMTO method, and excellent agreement 
with experimental data were obtained. fH^ ITH] In the following we summarise some details 
of the calculations. 

The dynamical matrix of fcc-Cs is calculated for a set of irreducible q points in a (8,8,8)- 
reciprocal lattice grid (29 points per l/48th part of the Brillouin Zone (BZ)). The (I,J,K) 
reciprocal-lattice grid is defined in the usual manner: qijk = {i/I)Gi + (j/ J)G2 + {k/ K)Gz, 
where Gj are the primitive translations in reciprocal space. We chose the exchange-correlation 
potential of Vosko-Wilk-Nusair|20j plus the non-local generalised-gradient-approximation 
(GGA-96) correct ion [21]. Because this gave a better prediction of the equilibrium volume 
at ambient pressure and OK which is only 4% smaller than the experimental zero-pressure 
volume at OK. On the other hand, the local density approximation (LDA) gave a nearly 20% 
overbinding at ambient pressure. The calculated bcc to fee transition pressure of 2.0GPa is 
in good agreement with the experimental value of 2.3GPa. jlj 

In the calculations a Sk — spd LMTO basis set is used with the one-center expansions 
inside the non-overlapping muffin-tin (MT) spheres performed up to Imax = 6. In the inter- 
stitial region, the s-, p- and d-basis functions are expanded in plane waves. The 5s semicore 
state is treated as valence state in a separate energy window. The induced charge densities 
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and screened potentials are represented inside the MT spheres by spherical harmonics up 
to Imax = 6 and by plane waves with a 148.4 Ry energy cutoff (9984 plane waves) in the 
interstitial region. 

The k-space integration is performed over a (16,16,16) grid (145 irreducible point) by 
means of the improved tetrahedron method^22 , but the integration weights for these k points 
are calculated from a denser (32,32,32) grid (897 points in the irreducible BZ). This results 
in a more accurate representation of the Fermi surface with a small number of k points. 

The phonons along the high-symmetry lines presented in the next section were calculated 
in a denser q mesh which fit to the (16,16,16) k-grid. 

3 Results and Discussions 

3.1 Phonon dispersions at V/Vq = 0.44 

Fig. 1 shows the calculated phonon dispersion curves along some high-symmetry directions 
for fcc-Cs at the volume 0.44Vo, which is in the experimental volume range of the cubic to 
tetragonal phase transition. To the right the calculated phonon density of states (DOS) is 
plotted. The dash lines represent the fitted results by the Born-von Karman model which 
is described in the next section. The calculated phonon frequencies at the high-symmetry 
zone-boundary (ZB) points L,X,K and W are listed in Table 1. 

As observed in Fig. 1, the phonon frequencies of the L[^^0], T[00^] and T[^^^] branches 
exhibit rather normal dispersion while the T[^^0], i^[00^] and i^fi^i^i^] dispersion curves show 
abnormal behaviour for intermediate and short wavelengths, which indicates long-range in- 
teractions between atoms. For the -^[^^^] phonon branch, small downwards anomaly is seen 
around ■C = 4- Near the zone boundary, the /^[■C^^] phonons have very little dispersion and 
the frequencies of the i^[00.^] phonons are rather small compared to other fee metals. The 
T[^^0] phonons have flat regions around ^ = | and | for the polarisation [llO] (Ti) and 
[001] (T2), respectively. Since the phonon DOS in Fig. 1 was calculated at a coarse (8,8,8) 
grid, these abnormal behaviours are not clearly seen in the plotted DOS. Nevertheless, the 
sudden increase of DOS at the frequency just above 0.3 THz and the small sharp peak at 
the frequency near 1.6 THz indicate the presence of these anomalies. 

In metals at equilibrium the long-range interatomic interactions are generally of little 
importance due to screening effect of the electrons. The lattice dynamical properties of met- 
als can therefore be described well using a short-range force model. This, however, may not 
hold for cesium under high pressure. The pressure-induced electronic s—>-d transition causes 
the valence charge density of cesium to deviate from spherical symmetry with a consequent 
reduction in electronic screening and increased interactions between the atoms. To check 
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this, we calculate the phonon dispersion curves of fcc-Cs at the volume K/Vq = 0.44 using a 
short-range force model within the quasiharmonic approximation. The dispersion relations 
are derived as in the usual harmonic approximation by considering non-central forces be- 
tween atoms up to second neighbours and the four interatomic force constants involved in 
the dynamical matrices are obtained from the three calculated elastic constants 23J and the 
T[100] zone boundary phonon frequency by expanding the elements of the dynamical matrix 
and finding its long- wavelength limit. The obtained dispersion curves are plotted in Fig. 1 
by solid lines and the phonon frequencies at the zone boundaries are listed in Table 1. Note 
that the quasiharmonic results presented here are derived from a short-range force model, 
which does not exactly consider the electronic contribution to the dynamical matrix. The 
difference to the results from linear-response calculations may therefore reveal the effects 
of long-range interatomic interactions. It may be seen in Fig. 1 that the short-range non- 
central force model describes the T[00^], T[^^^] and i^[^^0] phonon dispersions quite well. 
However, the previously mentioned anomalies in the L[00,^], L[^^^] and T[^^0] dispersion 
curves calculated by the linear-response method are not reproduced and they may therefore 
be a consequence of long-range interactions. 

In general, an abnormal behaviour of the phonon dispersion may be attributed to special 
properties of interatomic force constants. In the following subsection, we shall analyse the 
dispersion curves by using the Born-von Karman model to find these special properties of 
the interactions between the atoms in fcc-Cs. 

3.2 Born-von Karman analysis of force constants 

Since one of the motivations of the present work is to interpret the anomalies in the calculated 
phonon dispersion curves for fcc-Cs, it is necessary to use a model which can reproduce these 
anomalies. This could be achieved by a force- const ant model based on the Born-von Karman 
theory j211. 

In the Born-von Karman theory the frequencies i/(q) of the normal vibration modes for 
an fee lattice with one atom per unit cell are given as the solutions of a 3x3 determinantal 
equation 

|47r2Mz/2(q)<5„^-C,^(q)| = 0, (1) 

where M is the mass of the atom and 

CaM = E ^o.t3{l, l')exp[tii ■ R(/, /')], (2) 

I' 

where $a/3(/, /') is the ajS component of the force constant matrix between the atoms in the 

/th and I'th cell. For the phonon modes along the three high-symmetry directions [^^0], [00^] 
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and [^^^] , the C^p matrices are diagonal and the solution of Eq. ((T]) leads to equations in z/^ 
of the form 

A'k'^Mv'^ = ^ $„[! - cosin-nq/qmax)], (3) 

where qmax is half the distance to the nearest reciprocal lattice point in the direction of q, 
q = |q| and $„ is a sum of $«/?(/, /')ES] ^oi which the phase q ■ R(/, /') is a constant. $„ 
effectively represents a force between a plane of atoms and the planes of atoms normal to 
q and n planes away. The summation includes A^ terms so that $ri = for n > N. Thus, 
a Fourier series analysis of the squares of the phonon frequencies will yield the interplanar 
force constants. 

Accordingly, a least-squares Fourier fit are performed using Eq.© for the phonon dis- 
persion curves of fcc-Cs at the volume 0.44Vo calculated from linear-response calculations. 
The fitted dispersion curves are plotted in Fig. 1 by the dashed lines. Within the accuracy 
of the calculated phonon frequencies, it is found that a satisfied fit can only be obtained by 
including at least four planes in the [00^] branches, five planes in the T[^^^] branch and even 
six planes in the [^^0] and i^[^^^] branches. In Table 2 we list the fitted interplanar force 
constants $„. Since four planes along the [00^] direction correspond to interactions out to 
at least eighth neighbour atoms and five (six) planes along the [^C,^] ([^^0]) direction out to 
farther neighbours, the fitted results confirm the presence of long-range interactions between 
atoms in fcc-Cs. We plot the fitted $„ for the ^[110] [^^^] branch in Fig. 2a as a function of n, 
which effectively is the distance between planes of atoms. The $„ are observed to oscillate 
with the distance between planes, especially with a prominently negative $n=2- The values 
of z/^ for the same phonon branch and the fitted curve with A^ = 8 are shown in Fig. 2b 
together with the individual contributions from the first five Fourier components. It is seen 
that the n = 2 planes produce a largely negative contribution to the frequencies of phonons, 
which partly cancels the dominant contribution from the n = 1 planes. Consequently, the 
soft intermediate- wavelength T^iio] [^^0] phonons result and a flat region of the dispersion 
curve is formed. A similar oscillating behaviour of the fitted $„ is seen in some of the other 
branches although they are less pronounced. The anomalies observed in the T^ooi] [^^0] and 
-^[^'C'C] dispersion curves can be associated with these oscillations. For the i^[00.^] branch no 
oscillations in $„ is found. The very low phonon frequencies for long-wavelengths compared 
to the short-range model calculation may be attributed to the small value of $„=!• 

The interplanar force constants can be used to derive interatomic force constants |25j. 
With a general force model, a linear least-square fitting analysis for the $„ in Table 2 only 
allows us to work out the interatomic force constants in fcc-Cs extending to the fourth 
neighbour. The truncation of the interactions with farther neighbours introduces large er- 
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rors in the fitting procedure and makes the fitted interatomic force constants less rehable. 
Nevertheless, in Table 3 we list the best fit values of the interatomic constants in fee cesium. 

3.3 Soft mode and lattice instability 

The volume dependence of the high-symmetry zone-boundary L, X, K, and W phonons for 
fcc-Cs are presented in Fig. 3. From this the mode-Griineisen parameter ji = —{dlnvi/dlnV) 
was derived and listed in Table 1. It may be seen that all the high-symmetry phonons except 
the transverse X and L phonons decrease with increasing pressure and consequently have 
negative mode-Griineisen parameters, which may be an indication that the fee structure is 
unstable under pressure. 

In Fig. 4 we show the Tfiioji^^O] phonon-dispersion curves of fcc-Cs for the volumes 
V/Vq = 0.37,0.40,0.41, and 0.44. At volumes between V/Vq = 0.41 and 0.40 the phonon 
frequencies around ^ = | become imaginary. However, the shear elastic constant C does not 
become negative before the volume is reduced to V/Vq = 0.37 where the entire phonon branch 
becomes unstable. Christensen et al. fHl and Xie et al. flj inferred from their calculations 
that it is the pressure-induced negativity of C which is the cause of the instability of the 
fee phase. They however did not trace the development of the instability in details and it is 
clear from Fig. 4 that it is a qt^O T^iiq] [^^0] soft mode which drives the dynamical instability 
of the fee phase close to V/Vq = 0.41. As a consequence of the soft mode, superstructures 
should exist as a precursor to the tetragonal Cs-IV phase transition. It is interesting to note 
that the calculated transition volume V/Vq = 0.41 ~ 0.40 is close to the volume where the 
fee to tetragonal transition is observed [Sj. In Fig. 5 we illustrate the displacement pattern of 
atoms in fcc-Cs corresponding a Tjijo] [||0] phonon mode and the resulting Cs superstructure. 
As seen from Fig. 5, a structural element of the displacement pattern is triangular prisms, 
which are also present in the Cs-IV structure, even though other possible soft modes as well 
as the coupling between phonons could be important. 

Usually a soft phonon mode leads to a second-order or nearly first-order phase transi- 
tion. Early experimental results|5j indicated that the Cs-III-Cs-IV transition is first-order 
in character with a 4.3% volume reduction. Therefore soft phonon modes should not be 
considered the unique driving force in the cesium III-IV phase transformation even though 
the soft mode leads to an instability in the fee structure. The changes of the electronic 
properties under pressure, such as s — d transition, must play an important role. 



4 Conclusions 

Using a density functional linear-response LMTO method, we have studied the phonon spec- 
tra of fee cesium at volumes near the III-IV transition. Several anomalies in the dispersion 
curves have been found. 

Within the Born-von Karman theory of lattice dynamics, the observed anomalies in 
the calculated phonon dispersion curves for fcc-Cs at the volume 0.44Vo are connected with 
oscillating behaviours of the interplanar force constants against plane spacing n. In par- 
ticular, a negative (llO)-interplanar force-constant $„=2 is found to be responsible for soft 
intermediate-wavelength Tyiio] [^^0] phonons. 

The frequencies of high-symmetry K and W and longitudinal X and L phonons decrease 
with volume compression and the consequent negative mode-Griineisen parameters indicate 
an instability of fee lattice under pressure. When the volume is below 0.41Vo, the frequencies 
of the Tjiio] [^^0] phonons with ^ around | becomes imaginary prior to the sign change of 
the shear elastic constant C and the fee lattice is thus dynamically unstable. It is suggested 
that superstructures corresponding to the qT^O Tjijo] [■C'CO] soft phonon should be present 
as a precursor of tetragonal Cs-IV structure. It was pointed out in Ref. that the fee 
structure of Cs-III could not be confirmed and that Cs-III appears to be a rather complicated 
structure. The present argument may stimulate further experimental study for reexamining 
the structure of Cs at the pressure near the Cs-II-Cs-III-Cs-IV transitions. 
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Tables 



Table 1: Calculated phonon frequencies Vi (THz) at the high-symmetry L, X, K and W points for 
fee Cs at the volume V/Vq = 0.44 and corresponding mode-Griineisen parameters 7, defined by 
7j = —{dlnvi/dlnV). Also listed in brackets are the frequencies from the short-range force model 
calculation. 







L 


X 


K 


W 


L -branch 


^i 


1.97 (2.13) 


1.73 (2.10) 


1.58 (1.65) 


1.51 




li 


-0.95 


-0.79 


-0.80 


-1.27 


Ti-branch 


Vi 


0.82 (0.86) 


1.36 (1.36) 


0.95 (1.20) 


1.23 




li 


0.71 


0.49 


-0.80 


-1.60 


T2-branch 


Vi 


0.82 (0.86) 


1.36 (1.36) 


1.62 (1.93) 


1.46 




li 


0.71 


0.49 


-0.93 


-1.41 



Table 2: Interplanar force constants $„ (up to n = 6) for fee Cs at the volume V/Vq = 0.44 in 
units of 10"^ Ry/bohr^. The errors given are approximate errors for the fitted <&„. Note that 
a satisfactory fit for the [00.^] and T[^^^] phonons is obtained with only four and five Fourier 
components, respectively. 





'^i 


^2 


$3 


CI.4 


^5 


$6 


error 


r[ooe] 


10.69 


0.11 


-0.24 


-0.18 


— 


— 


±0.15 


L[m] 


16.09 


1.85 


0.72 


0.29 


- 


- 


±0.50 


Tm 


4.11 


0.49 


-0.36 


0.10 


-0.05 


— 


±0.02 


L[m] 


23.09 


-0.86 


-2.44 


3.35 


1.39 


-2.78 


±0.50 


r[ooi][ceo] 


16.32 


3.56 


-0.40 


-1.06 


0.80 


-0.55 


±0.19 


T[m][m 


8.51 


-4.10 


1.74 


-0.42 


0.21 


0.09 


±0.11 


L[m 


8.88 


14.52 


1.29 


-0.85 


0.30 


-0.29 


±0.12 
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Table 3: The best fit values of interatomic force-constants for fee Cs at the volume V/Vq = 0.44 in 
units of Ry/bohr^. The force constant notation follows the definition of Brockhouse et al.^2F. 



Neighbour location 


Force constants 


Values (Ry/bohr ) 


First 
§(1,1,0) 


ai 

71 



71 

ai 







/3i 


ai = 21.2x10-4 
/3i = -2.0x10-4 
71 =38.6x10-4 


Second 

§(2,0,0) 


02 







132 







/32 


a2 = -2.1x10-4 
/?2= 5.8x10-4 


Third 

§(2,1,1) 


as 

73 


S3 

P3 
73 


73 
73 
/33 


as= 3.6x10-4 
/J3 = -0.6x10-4 
7s = 3.6x10-4 
^3 = -0.4x10-4 


Fourth 
§(2,2,0) 


a4 

74 




74 

a4 






/?4 


04 = -0.4x10-4 
/34 = 0.04x10-4 
74 = -4.1x10-4 



11 



Figure Captions 

Figure 1: Calculated plio non-dispersion curves along three symmetry directions for fcc-Cs at 
the volume V jV^ = 0.44. The squares are the calculated results using the density-functional 
linear-response method and the solid lines represent the results from a quasiharmonic calcu- 
lation (see text). The dash lines are the fitted results by the Born- von Karman model. Also 
plotted in the right panel is the phonon density of states derived from the linear-response 
calculations. 

Figure 2: (a) The interplanar force constants (up ton = 8) for the T^iio] [CC^] branch phonons 
of fcc-Cs at the volume V/Vq = 0.44. A clear oscillating behaviour against the number of the 
plane (interplanar spacing) is seen, (b) Values of z/^ for the same branch phonons plotted as 
a function of the reduced wave- vector. The fitted curve (solid line) with 8 planes are shown 
together with the contributions from the first five Fourier components. 

Figure 3: Calculated phonon frequencies z/j (THz) at the high-symmetry L, X, K and W 
points for fee Cs as a function of the volume, (a) Longitudinal and (b) transverse branches. 

Figure 4: Calculated T^ijo] [■C'CO] phonon-dispersion curve for fee Cs at the volumes V/Vq = 
0.37,0.40,0.41, and 0.44. The points are the calculated values and the lines result from the 
interpolation between points and are only a guide to the eye. 

Figure 5: Schematic illustration of displacement pattern (left) of atoms in fcc-Cs correspond- 
ing a 7'[iio][4 40] phonon and the resulting Cs superstructure (right). The lengths of solid 
arrows are proportional to the displacement of atoms and the dash arrows indicate the [100], 
[010] and [110] directions of fee lattice. 
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